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Abstract 

We study the connection between charged lepton electric dipole moments, di (I = 
e, /i, r), and seesaw neutrino mass generation in a simple two Higgs doublet extension of 
the Standard Model plus three right-handed neutrinos (RHN) 2V , a = 1,2,3. For RHN 
with hierarchical masses and at least one with mass in the 10 TeV range we obtain the 
upper bounds of |cZ e | < 9 x 10 -30 e-cm and |cL| < 2 x 10~ 26 e-cm. Our scenario favors the 
normal mass hierarchy for the light neutrinos. We also calculated the cross section for 
e~e~ — > W~W~ in a high luminosity collider with constraints from neutrinoless double 
beta decay of nuclei included. Among the rare muon decay experiments we find that 
u — > e7 is most sensitive and the upper limit is < 8 x 10 -13 . 



1 Introduction 



In the celebrated seesaw mechanism pQ, the right handed Majorana neutrinos are essential 
to generating small Majorana masses for the active left-handed neutrinos of the Standard 
Model(SM). These fields are singlet under the SM gauge group and the exact number required 
is open to debate. Since the light neutrino masses are constrained to be less than an eV the 
masses of the right handed neutrinos have to be heavier then 10 12 GeV. This fits in well with 
expectations of grand unified theory although the seesaw scale is lower than the GUT scale 
which is generally taken to be around 10 16 GeV as required by proton stability. Moreover, 
the high scale also makes the seesaw mechanism impossible to test directly. The best we can 
hope for are indirect tests such as leptogenesis or renormalization effects. However, in order 
to make predictions in these latter studies additional assumptions have to be made and the 
results become highly model dependent. On the other hand, Majorana masses for the light 
active neutrinos can be tested in neutrinoless double beta decays of nuclei. Even in this case 
one has to eliminate other possible sources of lepton number violating new physics such as 
exotic scalars. Recently there are attempts to lower the masses of the right handed neutrinos 
to the TeV in leptogenesis studies Pj. They are particularly useful in supersymmetric models 
Clearly such low scale Majorana neutrinos are of phenomenological and theoretical interests 
in their own rights. They can be detected in high energy colliders and due to their rich CP 
properties effects in low energy experiments can also be searched for. The prominent example is 
the electric dipole moment (EDM) of a charged leptons denoted by di where I = e, /x, r. Already 
the experimental limit on \d e \ is an impressive < 10 _27 e-cm |3] and will be further improved in 
new round of experiments. In contrast, the limit on \d^\ < 10 -19 e-cm is much less stringent 
and dedicated experiments are now being proposed. 

In this paper we investigate the contribution of TeV scale Majorana right handed neutrinos, 
Nr, to d\. An immediate issue is to decide whether they are involved in generating active 
neutrino masses. Naively one expects that if they do so then the seesaw mechanism will 
restrict their Yukawa couplings to be very small and thereby making their contribution to d\ 
be minuscule. Thus, they can play an important role in the seesaw mechanism only in a subtle 
manner. Radiative effects on the left-handed part of the seesaw mass matrix was calculated 
in [3] . Previously it was pointed out [S] small Yukawa couplings can be avoided for more than 
one Nr. We shall display this in the context of a simple model which consist of the SM plus at 
least 3 right handed Majorana neutrinos, Nr, and an additional Higgs doublet. In order not to 
be confused with possible CP violating phases from the scalar potential we assume one Higgs 
doublet couples to Ir and the other to Nr. This is the natural flavor conserving extended two 
Higgs doublet studied in [7j and is also well known to be part of the minimal supersymmetric 
standard model. We shall see that one Nr can be arranged to be heavy and is responsible for 
the seesaw and the other two can be much lighter. Furthermore their Yukawa couplings can 
be of order unity. We do not attempt a detail fit to the neutrino mixing data which can be a 
separate study but merely to demonstrate the possibility of such a scenario. This very simple 
set up also gives rise to a nonvanishing di at the 2-loop level via a set of Feynman diagrams 
specific to Majorana fermions as pointed out by [HJ; a mechanism which has been checked by 
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|H] and JHl- The latter also contains a detail discussion of the 2-loop integrals. However, we 
will concentrate more on the structure of CP violation that appears and will be satisfied with 
order of magnitude estimates. We will be able to give a 'natural' order of magnitude estimate 
of the upper limit on d\ coming from the possible existence of multi-TeV scale Nr. 

One would wonder why do we add one more Higgs doublet in our construction. With only 
one Higgs doublet the 2-loop contribution to d\ from Majorana neutrinos is negligible [51 HUj. 
This is because now only the active neutrinos take part as they couple to the W bosons. Then 
di is proportional to their mass squared differences which are known to be small from neutrino 
oscillation data. Thus, in the SM extended to include seesaw neutrino mass although d\ happens 
at 2-loop it is still undetectably small. In contrast the SM with massless neutrinos d\ receives 
contribution at 4-loop or higher. With more than one Higgs doublet the physics changes. The 
right handed neutrinos do not decouple as we shall see later. This model can also serve as a 
prototype in studying the interplay between scalars and Majorana fermions in EDM's. 

In section 2 we describe in detail a model with three Nrs. We show how it works to generate 
sub-eV neutrino masses with one of them having mass in the 10 TeV range and he others can 
be much higher. For early pioneering work on Majorana neutrinos in gauge theories see [TT] . 

Next we discuss the Majorana phases in the model and dive into the d\ estimates. Since the 
physics scale we are interested in is relatively low the expected renormalization group running 
of the parameters are not very significant and we shall ignore them. Another possible test of 
this mechanism can be done using the reaction e~e~ — > W~W~; in the event that the right 
handed neutrinos are not kinematically accessible at LHC or the linear collider. This possibility 
has not been examined before and is discussed in section 4. No neutrino double beta decays 
(Ou/3/3) of nuclei are discussed here. Section 5 examines the possible tests in the rare decays 
H — > e7 and \i — > 3e using the constraint we found previously. Discussions of other low energy 
probes of Majorana phases can be found in [TJ]. Finally we give our conclusions. 

2 A Simple Model with Right Handed Neutrinos 

The model we study is the SM with two Higgs doublet denoted by 0i and 02 and 3 right handed 
neutrinos N a with a = 1, 2, 3. The terms in the total Lagrangian of interest to us are given by 



where 2 = 2O"20 2 an d indices i,j = 1,2,3. The hypercharge of Higgs doublets are Yi = 
Yi = 1/2. This is the simplest of two Higgs doublet models (2HDM). The details of the scalar 
potential V(4>i, 02) is not important for us and will not spell out here. Note that a Z 2 symmetry 
can be applied to fields so that the lagrangian is invariant under the transformation: 



+ 




(1) 



Li 



01 , 02 — > -02 , N 



a 



-N a , e R -> e R . 



(2) 
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In this example, the right-handed singlets only couple to 02 which is the simplest way to 
enforce natural flavor conservation. Other assignments to accommodate flavon models can also 
be applied. However, such details are not necessary for us. All the fermions are in the weak 
eigenbasis. M a b is a complex symmetric 3x3 mass matrix for the right-handed singlets and 
a, b = 1, 2, 3. Without loss of generality, M a b can be chosen to be diag{Mi, M 2 , M 3 } where the 
eigenvalues M a can be made real and positive. The Yukawa couplings and Q a are complex 
but the Higgs potential are taken to be real. In doing so the only source of possible CP violation 
come from the Yukawa couplings. 

It is worthwhile to examine the number of physical phases in this class of models. The 
discussion is clearest by first going into the charged lepton mass basis. In this basis we denote 
the Yukawa coupling by (' — > (. Now consider the general case with N right-handed singlets and 
n left-handed active neutrinos coupling as above. The neutrino mass matrix is a (n+N) x (n+N) 
matrix in a block matrix form: 



(3) 



v( 
v( T M 

where v( is a (n x N) Dirac mass matrix and v is a generic vacuum expectation value of the 
Higgs fields. For instance, for the Higgs doublet model of Eq.(JTJ) the following substitution 
should be made: v — > v sin/5 where tan (3 = Vijv\. Returning to Eq.Q, the lower block matrix, 
M, is the (N x N) Majorana mass matrix of the Nrs. It is symmetrical and complex and thus 
contain iV(iV + l)/2 phases. These can be completely absorbed by the U(N) mixing among 
the number N of Nr states. Another way to view this is to use the freedom of phase choice in 
Nr to rotate away the iV complex phases in the eigenvalues. Once that is done the phases of 
Nr are fixed. The phases of the active vl are still free. A phase redefinition of the vl will then 
remove n phases from (n x N) Yukawa terms. This leaves a total of n(N — 1) physical phases 
|13j . For the case of n = 3 and N = 3 we have 6 physical phases. It is customary to assign 
one phase to the light neutrino mixing matrix and leave the others in the mass eigenvalues. 
Moreover, the physics of EDM is seen more clearly using complex Yukawa couplings. 
The light neutrino Majorana mass matrix can be solved from Eq.(J3J): 



m u = -vXM'X ■ (4) 



2 / CilCjl . Ci2Cj2 Ci3Cj3 



Explicitly, the matrix elements are 

where i,j = e,fi,r. The standard seesaw mechanism is to assume Mi ~ M 2 ~ M 3 ~ 10 14 
GeV and the Yukawa couplings are all of order unity so as to get sub eV masses for the active 
neutrinos. It is also noted by many that it hard to obtain the observed bilarge mixing of the 
active neutrinos with inverse hierarchical masses. Without more assumptions we can extract 
one more result, i.e. 

AdetQ 2 , R x 
detmj = ————— . (6) 

1 1 M 1 M 2 M 3 v ; 
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which may be useful for constructing neutrino mass models. 

On closer examination of Eq.([5|) one discovers other ways of getting small neutrino masses. 
First we scale out the lowest of the three Nr masses which we call M < . Thus, Eq.flHJ) becomes 



where r a = M a /M <1 a = 1, 2, 3 and r a > 1 by construction. Each term in Eq. (f7jl can be view 
as a complex vector and it is a sum of three such vectors. If they form a triangle than the 
element vanishes. The smallness of the active neutrino masses can then be due to nearly closing 
of the complex triangle even with a value of M< in the TeV range. Similar techniques have 
been used to construct different hierarchies for the Nr to yield the experimentally acceptable 
mass matrices for m v P~2]. We leave aside the question of whether this is fining tuning or 
manifestation of approximate family symmetry of the heavy neutrinos. We take it to be purely 
phenomenologically motivated. 

As an example, we describe a scenario in which 3 Nr can generate sub-eV active neutrino 
mass but possesses the features that some of them are light enough, say ~ TeV, so that the 
seesaw mechanism can be amenable to testing in the near future. Assume that Q a obey the 
following relation: 



where 5 V is a small parameter we introduced. Then the seesaw mass contribution to the 
lights neutrino from the 3 Nrs nearly cancel among themselves. In the limit 5 U = 0, Eq.(jHJ) 
has a geometrical interpretation, i.e. the three terms viewed as vectors in the complex plane 
forms an equilateral triangle. Then a non-vanishing 5 U is a measure of the deviation from this 
configuration. In our example the £i side is (1 + 5 V ) longer than the sum of the other two 
and the small part left is responsible for small value of m v . As noted before the model has six 
physical phases and there are nine complex Yukawa couplings. For definiteness we will choose 
the three (us to be real. Clearly the deviation 5 U can be associated with either one of the M a 
and our choice is for simplicity of discussion. To obtain an acceptable mass matrix we need 
further assumptions. Taking a hint from the charged leptons we assume the couplings ( ea are 
such that ( ea <C ~ C™ — (a- Then a light neutrino mass matrix of the normal hierarchy 
type emerges 



It is well known that Eq. (JOJ) gives the observed bilarge mixing angles. Notice the scale of m v is 
set by Mi and the parameters 5 U and ( also play a crucial role in determining its magnitude. 
Moreover, it is important to note that the Yukawa couplings are complex as explicitly displayed 
in Eq.fjHJ). These are physical phases which will enter into EDM considerations. 

Interestingly we can extract more information about possible hierarchies in M a . There are 
three cases we can imagine: 




(7) 
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(a) (b) 




e u L (N a ) e e v L (N a ) e 

Figure 1: 1-loop diagrams which may give no-zero EDM. 

I. Mi ~ 10 12 GeV and it sets the scale for light active neutrinos. If we take 5 U ~ £1 ~ 0.1 
then this is sufficient to ensure sub-eV neutrinos. As seen in Eq.fjOJ M 2 and M3 play no 
role in determining m v \ thus they can be as light as a few TeV. However, Eq.fjHJ) dictates 
that (*L 2 will be very small and hence will not give a detectable di. 

II. Mi <C M 2 < M 3 with Mi ~ lOTeV. For light neutrinos masses in the sub-eV range Eq.fjUJ) 
requires the product 5 u (f < 10~ 9 . Superficially one would expect that 5 U ~ Ci ~ 10~ 3 . 
Since we do not have a theoretical basis for the values of these two parameters it is prudent 
to use experimental constraints. We shall see later that neutinoless double beta decays 
limit Ci < 0.1. From Eq.(|HJ) we see that even if M3 is of order 10 8 GeV the corresponding 
C will be of order unity. This is the interesting case for d\. Certainly we can lower M\ by 
simultaneously reducing 5 U or £i- 

III. Mi ~ M 2 ~ M 3 and they are in the TeV range. In this case, either all the Yukawa 
couplings are small or GIM-like cancellation with small mass splittings will have to take 
place. Either case no interesting di arises. 

The above discussion is sufficient to illustrate the connection between the seesaw mechanism 
and d[ we set forth to seek. Now we move on to the discussion of EDM. 

3 EDMs of Charged Leptons 

In models with charged scalars and Majorana neutrinos, d\ can begin at 1-loop. The relevant 
Feynman diagrams are given in FigO We shall employ the mass eigenstates in our discussions. 
The mixings between the light active neutrinos and the heavy iVVs are expected to be very 
small. Indeed from the simple case of one family seesaw this mixing is given by (v/M. For the 
parameter values discuss in case (II) above we estimate mixing between an active neutrino and 
the right-handed singlet which we generically called 9 to be < 10 -3 . 

As seen in Fig. (^1) both heavy right-handed neutrinos and light active neutrinos can enter. 
Moreover, the Yukawa couplings are conjugate of each other. Also the necessary helicity flip 
occurs in one of the external charged lepton lines; hence, there is no EDM from this diagram. 
For the W boson exchange diagram, see Fig^b), active neutrinos exchanges are dominant with 



5 



a small admixture of the right-handed neutrinos entering. In either case the two W vertices 
are also conjugate to each other and clearly there is no EDM from this diagram. 

At the two loop level there are 4 distinct topology we have to consider as shown in FigEl 
FigJ2fb-d) do not give rise to di when the thin lines represent gauge bosons. This is well 



(a) (b) (c) (d) 




Figure 2: The topology of 2 loop diagrams, the thick lines represent the fermion lines, and the 
thin ones could be scalars, gauge boson, or fermion loop (sub-diagram d). 

known form the SM. When they represent scalar particles the Yukawa couplings involve come 
in conjugate pairs ; thus negating their contributions. Since we have no phases in the scalar 
sector we arrive at the result that FigJ^b-d) give no EDM. This leaves only Fig|2fa) as the 
only type that can lead to a no n- vanishing d\. 

To see the physics more clearly we put the details in FigJ3J The external photon can attach 
to any charged object in the loops. The Majorana mass insertions for the light neutrinos are 




Figure 3: 2-loop diagrams which give non- vanishing EDM. 



indicated by the open box. These also flip helicity and change lepton numbers by 2 units. The 
corresponding insertions for the heavy Nrs are denoted by the filled boxes. The diagrams with 
the internal Uj lines replaced by iV a or vice versa are multiplied by the mixing 9 alluded in the 
discussions at the beginning of this section and will be suppressed. 

It is instructive to examine the two W-boson diagram. It has two open box insertions which 
involve light active neutrinos. Summations over different neutrino species and the three charged 
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leptons U are to be taken. The active neutrino mixing matrix elements at the incoming and 
outgoing lepton vertices can be different; thus leading to a non- vanishing di. The open box 
insertions indicating the lepton number violating nature of the Majorana masses are mandatory 
for this to happen. They do not exist for Dirac neutrinos and thus d\ cannot happen at the 
two loop level with Dirac neutrinos. Explicit calculations jHJ ITU] show that this diagram gives 
a contribution proportional to w? v and hence is completely negligible. Similarly, the diagrams 
FigJHfa) and (b) are suppressed by powers of m„/M w and also y e ~ 10~ 6 . These graphs can 
be neglected. This leaves only the two charged Higgs exchange diagram. The lepton EDM can 
be estimated as: 

di \ ^ m l ^ M a -M b M H 

~ ~ 2. (1^2)2 L /m ^^ (M B + M 6 )3 ln (M a + M b ) (10) 

a<b i=e,ti,T v ' v ' 

in the limit that the Nr are heavier that the charged Higgs boson which we assume to be of the 
weak scale. Strictly speaking, in the mass eigenbasis, the Yukawa couplings should be modified 
due to the mixing with the active neutrinos. These are expected to be small, i.e. 0(m u /M), 
and can be neglected. Note that the imaginary part of the product of four Yukawa couplings 
flips sign when one exchanges indices a <-> b. In other words, only the antisymmetric part in the 
loop integral yields the desired EDM operator. On the other hand, the factor (M a — Mb) also 
reflects that the CP violating effects go away when the masses of two right-handed neutrinos 
become degenerate. It is also interesting to note that the diagram with the photon attached 
to the internal charged lepton has no EDM contribution since the loop integral is completely 
symmetric in a and b. 

From Eq. fjlOj) we see that the EDM scales linearly as the mass of the charged lepton. For 
the hierarchical mass of case (II) and take M2 ~ M3 ~ 10 8 GeV and Mff ~ 200 GeV we obtain 
for the electron EDM 

\d e \ ~ 9.2 x lO-^lOTeV/MiflCei/O.lHCil 2 e-cm . (11) 

The above estimate is not very sensitive to the values of M# and since their dependence 
is logarithmic. Notice that we use a small value of C, e i as is required by the normal hierarchy 
solution and Ov/3/3 (see next section). If the right-handed neutrino masses are hierarchical such 
as Mi << M 2 < M 3 then we have 

\d e \ ~ 9.2 x 10- 31 (10TeV/Mi) 2 (M 2 /M 3 )|C e i/0.1| 2 |Ci| 2 e-cm (12) 

and will be suppressed compared to Eq. (fTT|) . Our estimate of d e is three orders of magnitude 
below current experimental limit JTo] but is within reach of new plan experiments ^E]- We 
note that Eq. (jlO|) is a good approximation to the actual two Feynman integrals which cannot 
be given in analytical form. Our numerical investigations show that it is accurate for order of 
magnitude estimates. 

We can also give an estimate of the muon EDM and it is 

|^| ~ 1.8 x 10- 26 (10TeV/Mi) 2 |Ci| 4 e-cm . (13) 
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In this case the internal r and /x lines give important contributions. When combined the 
coupling involves is d which need not be small. In contrast ( e i which is small enters the 
calculation for d e ( see Eq. (jl2j) ). Thus it is possible that can be more enhanced than 
just the mass factor m M /m e when compared to d e . The earlier discussion on the mass scaling 
violation of muon EDM is given in J7| . This illustrates the importance of doing both types of 
measurements. We add that our estimate is six to seven orders of magnitude lower than current 
limit PH] and will be a challenge even for the newly proposed dedicated d^ measurements ^H] ■ 

4 (0vf3f3) decay and its inverse 

Here we discuss how low scale Nr affects the decay rates of (Qvflfl) decays of nuclei. We will 
be concerned with the elementary quark level dd — > eeuu transition, and not worry about the 
detail nuclear physics. At the fundamental fermion level the amplitude is given by the diagrams 
below 

(a) (b) (c) 




Figure 4: Tree level (Ov/3/3) decay amplitudes from Majorana neutrinos 



An estimate of the amplitude for the 2W exchange diagram is 

A °~ 3 Mt~W <14) 

where (p) is the average momentum of the exchange light neutrino. The corresponding diagram 
with iV replacing v line is suppressed by the M and 6 2 . For the amplitude of Fig0Jb) with v 
exchange we estimate 

2 m g m e 1 m v,ee Q w i n— 11 A /ir\ 

4 ~»4iPIF~ !x " K ' (15) 



where m q represents the light quark mass, and with N exchange we obtain 

Ab ~ g 2 ™" ttStt-, }, ~ 3 x 10'^A a . (16) 

a A/f__. T\/r2 i\/r__ a v ' 
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Similarly the dominant 2H exchange graph is associated with a iV-line and it gives 



A, 
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2 m l 



M^M H M N 



6 x 10~ 8 A a 



;i7) 



In arriving the above estimation, we have used the following numerical numbers: ~ 10 
TeV, M w ~ lOOGeV, M H ~ 200 GeV, m Vfie ~ 10- 10 GeV, the quark mass m q ~ 10" 3 GeV and 
(p) ~ O.lGeV^ 1201 This analysis suggests that WW exchange with active Majorana neutrino 
exchange is still the dominated tree-level contribution to (0uj3j3) decay. 

Interestingly, for low scale right-handed Majorana neutrino important contribution to (Ov/3/3) 
decays can come from 1-loop diagrams depicted below 



e 




N , 




N ' 









W~ 



N 



N 




e~ H z 



W~ 



Figure 5: The box and triangle diagrams for (Oi//3/3) decay. 
The effective Lagrangian the diagrams generate is 

Fi (e c Lej e x ■ e 2 , 



£ 

Fi 



eg* 

32vr- 



E 



c 2 

M„ 



M 



H 



M„ 



where c is an order one constant which depends on the details of the 2HDM such as scalar 
mixings. The are polarization 4-vectors of the W bosons. We have also assumed M a 3> Mh 
with Mh denoting a common scalar mass. After dropping the m Ujee part the sum in the above 
equation yields 



a 



1 _2tH 

m^e 3 +lnr 3 e 



47TZ 

3 



(19) 



87rsin 2 ^Mi L" 

by using Eq.fjHJ) and assuming c = 1 for simplicity. If this is the dominant contribution to 
(Oz//3/3) decays then we can set the limit ( el < 0.1 for Mi = lOTeV. In comparison the low 
energy effective Lagrangian from the exchange of N± is given by 



C 



tree 



Ft (e c ~{^ u Le ) e£e% , 
Ana ( 6 2 



sin 2 6 



(20) 
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The numerical absolute value of the square bracket in Eq.(|19|) is ~ 5 for a large range of T2 
and r 3 . For £ el = .1 and M\ = lOTeV we get Fi = 6 x 10~ 9 whereas F t = 4 x 10~ n in units of 
GeV -1 . Thus, the loop diagram can be more important even when £ el is not that large. 

The effective Lagrangian of Eq. (jl8|) can also give rise to 2 W-boson production in e~e~ 
colliders even when N\ is too heavy to be directly produced. The cross section can be easily 
calculated to be 

a(ee _ ww) = JL_J_ (21) 

which is dominated by the longitudinal components of the W bosons. Although we have 
neglected the energy dependence in Fi this is sufficient for a ball park estimate of the cross 
section. For a linear collider with y/s = 2 TeV and a luminosity of 10 34 sec _1 cm~ 2 we obtain 1.3 
events in a 100 days running for a 10 TeV Majorana neutrinos. As mandated by the transient 
high s behavior of Eq. (j21j) one would require the highest available energy for a given high 
luminosity collider to probe this physics. It is easy to check that the usual tree level t-channel 
Nr exchange mechanism in the seesaw scenario gives a even smaller cross section 

5 A Trio of Rare Muon Decays 



V N a e 

Figure 6: Right-handed Majorana neutrino contribution to fi — > ej decay. 

The rare decays fi —>■ e + 7, 3e and fi — e conversion in nuclei have always been a favorite 
for testing models of lepton violations. If TeV scale Majorana iV/j's exist one also expects that 
these processes will occur. A general up to date review is given in (22 and we shall follow the 
notations used there. 

We begin with fi — > e'j. The most general Lorentz and gauge invariant fi — e — 7 interaction 
is given by 

M = -eAlW e { Pe ) { [f EQ (q 2 ) + f M o(q 2 h 5 ] lu (g Xu ' A/ 



q 2 



+ [fMi(q 2 ) + fEi(q 2 h 5 ] ^4 uM (22) 

where q x and A\ are the photon 4-momentum and polarization respectively and p e = — q. 
For fj, — > e7 only the form factors Jmi and Jei contribute. They can be calculated from the 
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dominant diagram given by Fig.©. The transition rate is 



fl f— ' >^%§r 8 - 01xl «-w(^) 4 - <*»> 

In arriving at the last formula, we have ignored the small Mu/M a term. The current experi- 
mental limit of < 1.2 x 10~ n j2B] sets a loose constrain on the mixings. Alternatively we can 
take Cei < 0.1 as required by our model and £ M i < 1 so as not to have strong Yukawa; then we 
get an upper limit of 8 x 10~ 13 for a 10 TeV Nr. We note in passing similar decays for the r 
are sensitive to the mixing £ rl which will be hard to obtain from other experiments. 

For /i — e conversion in nuclei the seesaw model belongs to the class where the photonic 
penguin diagram as given FigEl with the photon off shell dominates the transition rate. An 
explicit calculation gives 

ficonv = * / 9 §= e// Bfr - e + 7 ) • (24) 
IZir^L capt 

For f 2 Ti, F p = 0.55, Z eff = 17.61 and r capt = 1.71 x 10~ 18 GeV which implies B^ nv ~ 
0.0045 (// -> e + 7). 

Similarly for \x — > 3e the photonic penguin is the most important graph. The box diagram 
with two charged Higgs exchange is completely negligible. The Z-penguin graphs are also 
sub-dominant. Thus we obtain simply 

2a / m,, 11 \ , , . 

B - 3e = — In -H - — S // - e 7 25 
Sty \ m e 8 / 

or B(n — > 3e) ~ 0.006B(;U — > e7). Hence, /i — >■ e 7 and r — > /i(e)7 are the most important 
processes to probe TeV scale seesaw. 



6 Conclusions 

In this paper we have given a detailed study of the connection between seesaw neutrino mass 
generation and charged lepton EDMs. The 2HDM we employed is simple and it captures the 
physics clearly and succinctly. It is expected to be a crucial part of any elaborate embedding 
of the seesaw mechanism into a grand unified picture. As noted previously if all the right- 
handed neutrinos have very high masses, i.e. > 10 10 GeV then di will be undetectably small. 
We found that it is crucial to have at least one Nr have a mass in 10 TeV or slightly lower 
range. In addition, not all the Yukawa couplings can be suppressed as in the charged leptons. 
The charged lepton EDM arises from two loop diagrams involving Majorana neutrinos and the 
associated physical phases that have no counter parts in the SM with Dirac neutrinos. Under 
favorable choice of parameters we estimated that the upper limits are \d e \ < 9 x 10 _30 e-cm 
and \dfj\ < 1.8 x 10 -26 e-cm for a 10 TeV Majorana neutrino. The parameters involved are 
consistent with a normal mass hierarchy for the light active neutrinos. 
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Interestingly, a right-handed neutrino in the 10 TeV mass range may also be required for 
a successful leptogensis j2J- It is reasonable to expect that the mass of the lowest Majorana 
neutrino is in the 10 TeV range especially in the supersymmetry context. Moreover, the direct 
production of N\ is out of reach for high energy colliders under discussion. Even so we considered 
how a high luminosity e~e~ collider in the TeV range can still probe their existence via the same 
sign 2W~ production. Complementing this we calculated that the rare muon decay \i —>■ ey 
at the level of 10~ 13 is found to be sensitive test of the scenario we discussed. Similarly, one 
can contemplate the rare Z decays into r/i or re. The signatures are clean and unmistakable. 
At 1-loop this proceeds via a similar diagram as given in Fig.©. We estimate this to give a 
branching ratio of < 10~ 16 which is too small even for a Z-factory. 

Interestingly, the CP violation from the see-saw mechanism has negligible effect in the quark 
sector. This demonstrates clearly that the search of neutron and electron and muon EDM's are 
independent powerful probes of physics beyond the SM. 



7 Acknowledgement 

WFC wants to thank K. Cheung and C.Q. Geng for useful discussions. This work is supported 
in part by the Natural Science and Engineering Research Council of Canada. 



References 

[1] M. Gell-Mann, P. Ramond and R. Slansky, in Supergravity ed. by P. van Nieuwenhuizen 
and D. Freedman (Noth-Holland 1979) 315; 

S.L. Glashow, in Quarks and Leptons Cargese Lectures, ed. M. Levy et al, (Plenum, New 
York 1980) 707; 

T. Yanagida in Proc. of the Workshop on the Unified Theory and the Baryon Number in 
the Universe ed. O. Sawada and A. Sugamoto (KEK Repeort 79-18) 95; 
R.N. Mohapatra and G. Senjanovic, Phys. Rev. Lett 44(1980) 912. 

[2] M. Fukugita and T. Yanagita, Phys. Lett., B174(1986) 45; 
A. Pilaftsis, Phys. Rev. D 56, 5431 (1997). 

[3] L. Boubekeur, T. Hambye, and G. Senjanovic, Phys. Rev. Lett. 93 (2004) 111601; 
T. Hambye, J. March-Russell and S.M. West, JHEP 0407 (2004) 070; 
J. March-Russell and S.M. West, Phys. Lett. 59 3 (2004) 181; 
M. Raidal, A. Strumia and K . Turzyhski, ||hep-ph/04080l"5l 
S.M. West, hep-ph/0408318, 



[4] M.V. Romalis, W.C. Griffith and E.N. Fortson, Phys. Rev. Lett. 86(2001) 2505. 
[5] A. Pilaftsis, Z. Phys. C 55, 275 (1992). 

[6] L.N. Chang, D. Ng and J.N. Ng, Phys. Rev. D50(1994) 4585; 

W. Loinaz, N. Okamura, S. Rayyan, T. Takeeuchi and L.C.R. Wijewardhana, idid 
D68(2003) 073001. 



12 



[7] S.L. Glashow and S. Weinberg, Phys. Rev. D15 (1977) 1958. 

[8] D. Ng and J.N. Ng, Mod. Phys. Lett., A 11 (1996) 211. 

[9] A. Pilaftsis, Int. J. Mod. Phys. A 14 (1999) 1811; 

A. Pilaftsis and T. E. J. Underwood, Nucl. Phys. B 692, 303 (2004). 

[10] J. P. Archaumbault, A. Czarnecki and M. Pospelov, |hep- ph / 0406089 , 

[11] J. Schechter and J.W.F. Valle, Phys. Rev. D 22 (1980) 2227; ibid D 23 (1981) 1666, and 
ibid D 24 (1982) 774. 

[12] A. de Gouvea, B. Kayser and R.N. Mohapatra, Phys. Rev. D67 (2003) 053004. 

[13] A. Santamria, Phys. Lett. B305 (1993) 90. 

[14] S. Antush and S.F. King, New J. Phys. 6 (2004) 110. 

[15] B.C. Regan, E.D. Commins, C.J. Schmidt, and D. DeMille, Phys. Rev. Lett. 88 (2002) 
071805. 

[16] M.G. Kozlov and D. DeMille, Phys. Rev. Lett. 89 (2002) 133001. 

[17] K. S. Babu, B. Dutta and R. N. Mohapatra, Phys. Rev. Lett. 85, 5064 (2000); 
J. L. Feng, K. T. Matchev and Y. Shadmi, Nucl. Phys. B 613, 366 (2001). 

[18] R.McNabb, hep-ex/04070"08| 

[19] F.J.M. Farley et al, Phys. Rev. Lett. 93 (2004) 052001. 

[20] H. Pas, M. Hirsh, H.V. Klapdor-Kleingrothaus, and S.G. Kovalenko, Phys. Lett. B453 
(1999) 194; 

J. Vergados, Phys. Rep C361 (2002) 1. 

[21] T.G. Rizzo, Phys. Lett. B116 (1982) 23; 

D. London, G. Belanger, and J.N. Ng, Phys. Lett. B188 (1987) 155. 

[22] Y. Kuno and Y. Okada, Rev. Mod. Phys. 73 (2001) 151. 

[23] Particle Data Group, Phys. Lett. B592 (2004) 409. 



13 



